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The representation of a quantum system as the spatial configuration of its con¬ 
stituents evolving in time as a trajectory under the action of the wave-function, is the 
main objective of the De Broglic-Bohm theory (or pilot wave theory). However, its 
standard formulation is referred to the statistical ensemble of its possible trajectories. 
The statistical ensemble is introduced in order to establish the exact correspondence 
(the Born’s rule) between the probability density on the spatial configurations and 
the quantum distribution, that is the squared modulus of the wave-function. In this 
work we explore the possibility of using the pilot wave theory at the level of a single 
Bohm’s trajectory, that is a single realization of the time dependent configuration 
which should be representative of a single realization of the quantum system. The 
pilot wave theory allows a formally self-consistent representation of quantum sys¬ 
tems as a single Bohm’s trajectory, but in this case there is no room for the Born’s 
rule at least in its standard form. We will show that a correspondence exists be¬ 
tween the statistical distribution of configurations along the single Bohm’s trajectory 
and the quantum distribution for a subsystem interacting with the environment in a 
multicomponent system. To this aim, we present the numerical results of the single 
Bohm’s trajectory description of the model system of six confined rotors with random 
interactions. We find a rather close correspondence between the coordinate distribu¬ 
tion of one rotor, the others representing the environment, along its trajectory and 
the time averaged marginal quantum distribution for the same rotor. This might 
be considered as the counterpart of the standard Born’s rule when the pilot wave 
theory is applied at the level of single Bohm’s trajectory. Furthermore a strongly 
fluctuating behavior with a fast loss of correlation is found for the evolution of each 
rotor coordinate. This suggests that a Markov process might well approximate the 
evolution of the Bohm’s coordinate of a single rotor (the subsystem) and, under this 
condition, it is shown that the correspondence between coordinate distribution and 
quantum distribution of the rotor is exactly verified. 
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I. INTRODUCTION 


Since its origin, quantum mechanics gave rise to vivid debates about the possible physical 
interpretations of its mathematical formalism. The widely accepted Copenhagen interpre¬ 
tation 1 ^ associates the wave-function to the probability of the outcomes in the measure of 
observables of the quantum system by invoking a purely classical measurement apparatus. 
It represents an epistemological approach since it does not attribute physical meaning to 
the wave-function independently on the act of measure. In this framework the Born’ rule 
identifies the probability density for a particular configuration of the constituents with the 
square modulus of the wave function. 

In 1952 David Bohm elaborated a different interpretation 2 ® called “pilot wave theory”, 
previously suggested by De Broglie at the Solvay congress in the 1927 4 . Indeed it is called also 
“De Broglic-Bohm theory”. The pilot wave theory assumes that the physical constituents 
of the quantum systems are particles possessing like in classical mechanics well defined 
coordinates, i.e., the configuration of the system. The wave function plays the role of a 
held, like the electromagnetic held but in the configuration space, that pilots the particle’s 
motion. Thus, the particle’s coordinates and the wave-function are supposed to be reality 
elements, in other words an ontological status is attributed to both of them. The particle’s 
configuration together with the wave function represents the dynamical state of the system 
which evolves in time deterministically. In such a framework the system’s physical properties 
are specihed according to the particle’s configuration without reference to the measurement. 

This description in terms of the time dependent wave-function and of the trajectory of 
particle’s configuration together, is well defined formally but it leave open the issue whether 
a relation exists with other interpretations of quantum mechanics. In order to provide an 
answer, Bohm moved from a single trajectory picture of the reality to a representation with 
an ensemble of particle’s configurations, i.e., a swarm of trajectories evolving in time under 
the action of the same pilot, that is the same wave-function. A density distribution is 
defined for such an ensemble to describe the probability density that the set of particles is in 
a given configuration. Bohm introduced the postulate that the distribution of the system’s 
configuration at the initial time coincides with the square modulus of the wave-function 2 . 
Then one derives the equivalence at all times between the distribution of configurations and 
the square modulus of the wave-function, this being equivalent to quantum hydrodynamics 
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formulation of the Schroedinger equation by Madelung-. In this way, the agreement with 
the Born’s rule of standard quantum mechanics is assured. 

On the other hand there are no evident reasons for supporting on a physical ground the 
equivalence between the square modulus of the wave-function and the density of system’s 
configurations at the initial time. Bohrn himself recognized the critical role of this assump¬ 
tion and argued that if this is not the case, then the randomness deriving from particle’s 
interactions would enforce such a correspondence during the time evolution 6 . Such a point 
of view has been further developed recently by Valentini and coworkers with the objective 
of demonstrating that an initial arbitrary distribution on the configuration space relaxes in 
the time to the square modulus of the wave-function 7 9 . A different procedure has been 
proposed by Durr and coworkers by introducing an effective wave function representative of 
a system interacting with the environment 10 . 

The De Broglie-Bohm theory found applications in different aspects of quantum physics: 
the most common examples concern the analysis of specific phenomenas, as tunneling or 
scattering 11 ^ modeled by a wave packet motion of one particle systems. These studies 
show that a simple picture of particle’s trajectories can be derived for the quantum dy¬ 
namics. In other research fields the pilot wave theory has mainly a computational role in 
order to reconstruct the time dependent wave-function from a collection of evolving system’s 
configurations 13 15 . Furthermore pilot wave theory has been considered as the appropriate 
framework to address the semi-classical approximation of quantum mechanics 161 '. There 
are also attempts to use it as tool in the development of new multiscale procedures for 
large size systems, with traditional quantum computational approaches handling only a re¬ 
stricted member of degrees of freedom, while the remaining are treated according to classical 
formalisms 1 ®. 

In this work we intend to explore the application of the pilot wave theory at the level of a 
single trajectory of the system’s configuration. The underlying motivations derive from the 
recent investigations on single-molecule or single-spin observables 19 20 . This together with 
the efforts towards the realization of quantum computers based on nanostructures®®, calls 
for a representation of material systems according to a single realization of the quantum 
state. In this respect, the formal structure of the De Broglie-Bohm theory is well defined 
and self-consistent also when an unique wave-function and an unique time-dependent sys¬ 
tem’s configuration are used to describe a particular realization of the quantum system. 
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This point of view might be considered as the most natural way of interpreting the pilot 
wave theory without requiring any particular constraint on the distribution of initial config¬ 
urations. Indeed, the original formulation of the Bohm theory includes implicitly the idea 
that a real system has always a well defined configuration, and we intend to explore the 
implications of this picture of the real systems through a single trajectory. 

In this way, however, a major issue remains open: what is the connection with the 
Born’s rule of standard Quantum Mechanics? More specifically, how to define from a single 
trajectory of the system’s configuration a probability density on the particle’s coordinates, 
which is a prerequisite before to establish a relation with the wave-function? One can exploit 
the analogy with classical statistical mechanics which introduces the equilibrium distribution 
by considering the density of phase space points along a single time dependent realization of 
the isolated systemP 4 . Also in the case of the pilot wave theory with a single Bohm trajectory 
(i.e., the time dependent system’s configuration), one can define in an analogous way the 
equilibrium probability density on the set of particle’s coordinates. But then, how to compare 
such an equilibrium distribution with the quantum distribution given as the square modulus 
of the wave function, which is intrinsically a time-dependent distribution? As a matter of 
fact the comparison becomes meaningful when the marginal distributions are considered 
for a subsystem interacting with a larger environment acting as thermal bath. Indeed, by 
employing the methods developed in a previous work 2 ^, one can show that in such a case 
the fluctuations of the marginal quantum distribution become negligible. Our conjecture is 
that in this particular situation the marginal distributions obtained from the configuration 
distribution along a Bohm trajectory and from the wave-function tend to coincide. In order 
to provide evidences about this behavior, we shall present some computational results for a 
model system of several, randomly coupled, confined rotors. It should be evident that such 
a conjecture, if verified, plays the same role of the Born’s rule in the Bohm analysis of the 
ensemble of trajectories. 

Few attempts has been done to connect the Born’s rule with a single Bohrn’s trajectory. 
Shtanov 26 investigated the problem from the point of view of ergodicity. Very recently, 
Philbin 2 -^ considered a simple one dimensional system (an harmonic oscillator) in the pres¬ 
ence of an external time dependent potential which mimics the position measurement. From 
a temporal sequence of these position measurements he obtains the same distribution given 
by the square modulus of the wave-function. In spite of the differences on the employed 
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model systems and on the type of dynamical regime, we share the same objective of devel¬ 
oping the pilot wave theory for a single realization of the quantum system. 


The paper is organized as follows. Since we consider a model system composed of several 
interacting components, statistical tools are required to analyze the quantum pure state 
represented by the wave-function. In the next section we introduce the Random Pure State 
Ensemble (RPSE) employed for the sampling of the wave-function, and we summarize its 
fundamental properties 25 2830 . Such a statistical ensemble allows one to evaluate the ampli¬ 
tude of fluctuations of the quantum observables (expectation values) with respect to their 
equilibrium values defined by time averages, and to estimate the behavior of fluctuation 
amplitudes in the thermodynamic limit for increasing size of the system. By recalling the 
results reported in Ref. 2 ^, it is shown that the marginal distribution on a subsystem, as 
obtained from integration of the square modulus of the wave-function on the environment 
degrees of freedom, is characterized by fluctuations of vanishing amplitude for increasing size 
of the environment. Therefore the subsystem is described by a nearly stationary marginal 
quantum distribution if the environment is large enough, and it can be well approximated 
by its time average. In the following Section |III| the standard form of the pilot wave theory 
is summarized, and the procedure for generating the single Bohm’s trajectory is illustrated. 


In Section IV the model system of six interacting confined rotors is used to verify our conjec¬ 
ture. First the model system is described in detail together with the numerical procedures 
employed for the calculation of time dependent properties. Then the main results are illus¬ 
trated in relation to: 1) the nearly stationarity of the marginal quantum distribution of one 
rotor, the other five rotors constituting the environment, 2) the randomly fluctuating behav¬ 
ior displayed by the evolution of the Bohm’s coordinate of one rotor with the corresponding 
loss of correlation with the time, 3) the close correspondence between the marginal quantum 
distribution and the distribution of Bohm’s coordinate of the subsystem, which provides the 
computational evidence of our conjecture. In Section [V] we show that the conjecture is ex¬ 
actly verified if the Bohrn coordinate of the subsystem behaves like an independent Markov 


stochastic variable, as partially suggested by the numerical results. In the final Section VI 
the general conclusions are drawn by focusing on the implications of our work. 


6 






II. STATISTICS OF QUANTUM PURE STATES 


In this section we present the statistical description of quantum pure states to be employed 
in the analysis of single Bohm trajectories. The need of a statistics of pure states, that is 
of quantum states described by a wave function belonging to the proper Hilbert space for 
an isolated (closed) system, arose mainly from the efforts of demonstrating the typicality 
of quantum observabled 31 32 . On the other hand, well defined statistical rules are required 
for sampling the initial wave function whenever the quantum dynamics is examined without 
particular a priori choices of the initial state. We stress that in quantum mechanics the 
condition of isolated system is more stringent than in classical mechanics: entanglement 
would keep the system connected to his environment even though there is no energy exchange 
between them. 

In standard quantum mechanics 3 ® the wave function |\P(£)) is an Hilbet 7~L space’s element 
representing the state of the isolated system, which evolves in time through the Schrodinger 
equation from a given initial state |T(0)), 

= e~ lflt/h \y(0)) ( 1 ) 

H being the Hamiltonian of the system. The time evolution A[t) of a generic physical 
property described by self-adjoint operator A is determined by the expectation value, 

m = <«'W|i|S'(i)> = r rv{M*)}, (2) 

where g(t) is the density matrix operator for the pure state 


e{t) = \m)m) i- 


( 3 ) 


An expectation value A{t) is usually interpreted as the mean value of a infinite number of 
measures of the observable at time t. Like in classical statistical mechanics 2 ®, we identify 
the equilibrium value of observable A with the asymptotic time average of the expectation 
value: 


A(t) := lim — 

V ' T— >+ oo T 


d t A(t ) = Tr 


( 4 ) 


where g(t) is the time average of the density matrix. 

If we are interested to properties of a subsystem, we can imagine to partition the isolated 
system: the subsystem S, also denoted as system in the following, and the environment E 
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for the remaining part of the isolated system. Correspondingly the overall Hilbert space H 
is factorized into the Hilbert spaces Ws of the system and He of the environment, H = 
Hs®He- In this situation the observable a{t) of interest, i.e., a property of the subsystem, 
is represented by the expectation value of an operator a <g) 1®. The reduced density matrix 
operator <f(t), as obtained by partial trace Ti'e over the environment states of the pure state 
density operator, 

a(t) := Tr E {£(£)}, (5) 

allows the calculation of this expectation value within the subsystem Hilbert space 

a(t) = Tr|(a <g) i E )^(t)| = Tr 5 |a d(f) j. (6) 

Its equilibrium value a{t) is defined again by time averaging and it can be evaluated like in 
Eq. (Jr]) by means of the time average a(t) of the reduced density matrix. 

In order to formulate a statistical description of quantum pure states, a finite set of 
parameters identifying the instantaneous wave function has to be selected, very much like 
for the phase space of classical statistical mechanics. This requires the confinement of the 
wave function to a finite dimensional subspace of H, say a A r -dimensional subspace Hn in 
the following called as active space. To select the active space, it is convenient to resort to 
the orthonormal eigenstates | E/.) of the Hamiltonian: 

H \E k ) = E k \E k ) . (7) 

Like in previous works 28 ®Ell, we shall employ the following type of active space 

Un := {© \ Ek ) with En E rnax < E_\r ] J, (8) 

that is the subspace due to eigenstates with eigenvalues smaller than E max . The energy 
cutoff E max is the only parameter required for the identification of this active space, and 
it has been shown that in the limit of macroscopic systems E max represents the internal 

energ;y®. 

It should be mentioned that one can employ an alternative active space by introducing 
also a low energy cutoff E m j n , like in the definition of microcanonical density matrix of 
standard quantum statistical mechanic^ 5 . In this way, however, one has to manage two 
different cutoff parameters and, furthermore, no direct relation exists between the lower 
energy cutoff E min and thermodynamic propertied. 






The wave function | *!/(£)) at a given time is then conveniently specified as a linear 
combination of the basis elements \E k ) of the active space through its expansion coeffi¬ 
cients c k {t) or, equivalently, through the sets of populations (Pi,P 2 , ..., P/v) and of phases 
a(t) = (aq(t), ca 2 (t), ..., ceN(t)) obtained from the polar representation of the expansion 
coefficients 

c k (t) := (E k \*(t)) = s/p k e-^\ (9) 

with a linear time dependence of the phases: a k (t) = ccfc(O) + E k t/h. Because of the 
normalization condition, 

N 

<*wi* , w> = E p * = 1 ’ < 10 ) 

k= 1 

only (TV — 1) populations are independent, say the set P = (P 1; P 2 ,..., P/v_i). Therefore a 
bijection exists between the normalized wave function and the ensemble of populations P 
and of phases a(t), with each particular set of these (2TV — 1) real parameters corresponding 
to a specific wave function |\P(i)). In other words, all the pure states of the active space 
can be imagined like unit vectors drawing an unit sphere in a 2TV-dimensional Euclidean 
space®. 

Because of the choice of expanding the wave function along the Hamiltonian eigenstates, 
the phases are the only dynamic variables of the system, while the populations represent the 
constants of motion. Correspondingly it is easily shown that equilibrium properties A(t) of 
Eq. ([4]) depend on populations only. Indeed, under the condition of rational independence of 
Hamiltonian eigenvalues, meaning that equation Ylk=i n k E k = 0 for integer n k has only the 
trivial solution n k = 0 VP, the equilibrium density matrix is diagonal with the populations 
as diagonal element^®, 

N 

g(t) = ^2 Pk '> ( n ) 

k =1 

— p 

where we have introduced the symbol g to highlight the dependence of equilibrium density 
matrix on populations only. We emphasize that the condition of rational independence is 
not too restrictive because of the contribution of random interactions, typical of material 
systems, leading to energy eigenvalue distribution having, at least partially, a random char¬ 
acter 36 . Therefore, according to Eq. Q, also the equilibrium value of a generic observable 
depends on populations only, 

W) = Tr{i| P } = H P , (12) 
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—P 

as stressed by the symbol A . 

In this way, the parametric dependence of the equilibrium properties on the populations 
results to be evident. On the other hand, there are no empirical methods leading to a com¬ 
plete characterization of the initial state |T(0)) and, therefore, also of the populations. This 
implies that populations can be characterized only on a statistical ground by selecting the 
ensemble for their probability distribution. The absence of privileged directions for |\h(0)) 
within the unit sphere, leads quite naturally to a purely random choice for the ensemble of 
pure states. In Ref. 28 the Random Pure State Ensemble (RPSE) for populations has been 
characterized from the geometrical analysis of the measure on the unit sphere, so deriving 
the probability density on the (N — 1) independent populations P 

Prpse(-P) = (N — 1)!. (13) 

Such a probability density allows the explicit calculation of the ensemble average of equilib¬ 
rium properties, 

:= j dP x ... dP/v-i ~A P Prpse(Pi, • • •, Pn- i), (14) 

— P 

which can be interpreted as the average of A amongst random realizations of the initial 

pure state |T(0)). Notice that integration domain on populations is bounded by constraints 

0 < P fc < 1, \/k = 1,2,... N. In the following we shall employ the bracket ^ ^ to denote 

the RPSE average of a function of populations. 

In order to recover also the macroscopic description of the system, one should consider 

the equilibrium energy, H P = 'p!^ =1 PkEk, and Shannon’s entropj^ with respect to the 

populations, S p = —kp JJkLi Pk l' 1 Pk- Their RPSE average are associated respectively to 

- P 

the thermodynamical internal energy, U := (H ), and to the thermodynamical entropy, 
S := ( S p ), both being functions of E max . By eliminating the E max dependence between 
functions U(E max ) and S(E max ), one recovers the thermodynamical state function S(U) 
and the temperature as well from its derivative 1/T = dS/dU. In this framework, by 
considering the system as an ensemble of n distinct components, like molecules in material 
systems, one can define the thermodynamic limit for n —>• oo at a given temperature 30 . The 
thermodynamic limit requires the tensorial product of the Hilbert spaces of all the distinct 
components, and this implies an exponential growth of the dimension N of the active space 
"Hat with the number n of components 25 . Finally, in the same limit, the RPSE average of the 
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equilibrium reduced density matrix y a y of a subsystem having weak enough interactions 
with the environment, takes the canonical form 



e ~H s /k B T 

Tr s {e-'W‘s T } ’ 


(15) 


where Hs is the Hamiltonian of the subsystem 30 . 

The RPSE statistics allows the quantitative analysis of typicalitji 30 of an equilibrium 

—p 

property A by evaluating the thermodynamic limit of its square variance within the en¬ 
semble, 

^ P ~<^>) 2 )' <“> 

—P 

Typicality of property A is assured if this limit vanishes, this implying that the value of 

— P 

A in a realization of the pure state is independent of the set of populations, as long as its 
deviation from the ensemble average ^ A P ^ tends to vanish. In other words, property A P is 
typical in the meaning that it is nearly independent of the particular realization of the pure 
state. 

Furthermore, RPSE ensemble allows the quantitative analysis not only of typicality of 

an observable, but also of its time fluctuations which are of primary importance for the 

objectives of this work. In order to characterize the amplitude of fluctuations of A{t) during 

its time evolution, we consider the equilibrium value, i.e., the time average, of the square of 

- P 

deviations A A(t) := A(t) — A from the time average 

(KW P (A(t)-A F ) 2 (17) 


that, like all the equilibrium properties, depends on the population set. The population 
average within RPSE provides an estimate ^(AA) 2 ^ of squared fluctuations which is in¬ 
dependent of the particular realization of the pure stat^ 5 38 and reads 

<^ P ) + (( aP -( 3P ))) = FT7' < 18 > 

where the second term at the left hand side describes the typicality of equilibrium property 

- P 

A as previously discussed. At the right hand side, N is the dimension of the active space 
"Hat, while D 2 (A) represents the squared spectral variance of the operator A, D 2 (A ) = 
J2k=i(^k ~ Di(A)) 2 , where {A*,} is the ensemble of eigenvalues of A in l-i N and Di(A) = 
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Ylk=i is the eigenvalue average. Such a relation connects the statistical properties of 

the expectation value A(t), at the left hand side of the equation, to the spectral properties of 
the operator A, on the right hand side of equation. If operator A has a bounded spectrum, 
then D 2 (A ) is finite and in the thermodynamic limit, n —y +oo, the right hand side of 


Eq. (18) vanishes because of the exponential growth with n of the active space dimension 


N. Correspondingly also both terms at the left hand side of Eq. (18) vanish since they are 
non negative 


lim 

n—>•+oo 


A P -(A P 


= lim 

n —>-+oo 


(A^R) = 0. 


(19) 


Thus, in the thermodynamic limit, both typicality and the vanishing of fluctuations are 
assured for bounded operators. Outside the thermodynamic limit, for finite but large enough 

_ p 

isolated quantum systems a nearly stationarity A(t) ~ A is predicted. Furthermore, we 
note that in this conditions the expectation value A(t) is nearly equal to the thermodynamic 
value ^ A P ^ because of typicality: A{t) ~ ^ A P ^. 

These results for typicality and fluctuation amplitude of bounded operators can be applied 
to the reduced density matrix of a subsystem of an isolated system. In particular, as shown in 
detail in Ref. 2 ^, the following condition for the expectation value a(t) of subsystem operator 


a derives from Eq. (18), 


a p - (a p )^j ^ + <^(Aa)2 J 


< 


Trs< a 2 (d 


— p 


Tr 5 |a^o- P ^) 


( 20 ) 


N+l 

In the thermodynamic limit the ensemble average of the reduced density matrix tends to 


the canonical form Eq. (15) and, therefore, the right hand side vanishes because of the 


active space dimension N of at the denominator. Then both typicality and the vanishing 


of fluctuations are recovered like in Eq. (18) for bounded operators, but now for a generic 


operator a of the subsystem. For finite but large enough isolated systems this implies that 
subsystem observables are nearly stationary, 


a(t) ~ a 


-p 


( 21 ) 


that is, their time dependent deviations from the equilibrium values is negligible. 

As an application of the previous analysis, we examine the statistical distribution on the 
coordinates q$ for the subsystem degrees of freedom. In standard quantum mechanics the 
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wave function allows the calculation of the time dependent distribution on the ensemble of 
coordinates q = {qk}k=i,...,n of the isolated system with n degrees of freedom through the 
probability density 

p(Q,t) = \^(q,t)\ 2 (22) 

with a parametric dependence on the initial pure state determining the time dependent wave 
function. Once the subsystem S, and the environment E as well, has been selected, the 
isolated system coordinates can be identified with the ensemble q = ( qs,qE ) of subsystem 
coordinates qs and of coordinates qE for the environment degrees of freedom. Then, by 
integration on the environment coordinates, the marginal distribution on the subsystem 
degrees of freedom is recovered 

P S (qs,t ) := J d q E p(q S ,qE,t). (23) 

Like for any time dependent observable, the time average defines the corresponding equilib¬ 
rium property, in this case the equilibrium distribution 

p S ’ eq (q s ) :=p s (qs,t), (24) 

where the reference to the parametric dependence on population set P has been omitted for 
the sake of a compact notation. 

Let us consider now an orthonormal basis for the subsystem Hilbert space "Hs, 

and its representation {p m (qs)} as explicit functions of subsystem coordinates qs . For any 
set of qs values, we can define the following operator 

a(qs) ■= l Pm( ( ls)<Pm'(qs) » (25) 

m,m' 

where its operator nature is determined by the bras (tp m '\ and kets on the r.h.s.. One 
can easily verify that its gs-dependent expectation value supplies the subsystem marginal 
probability density calculated at qs 

P S (qs,t) = Tr s |% s )6-(t)}. (26) 

In this way, the marginal distribution can be interpreted as expectation value of a subsystem 
operator, which is characterized by typicality and absence of fluctuations in the thermody¬ 
namic limit in agreement with the previous conclusions. Outside the thermodynamic limit, 
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but for large enough isolated systems, negligible contributions of fluctuations about the time 


average Eq. (24) are expected like in Eq. (21), 


p s (qs,t) -p s,eq (qs), 


(27) 


so that the subsystem is characterized by a nearly time independent marginal distribution. 

In conclusion, as long as expectation values or, equivalently, marginal distributions de¬ 
rived from the wave function are employed to describe a subsystem which is part of a much 
larger isolated system, the time evolution of the subsystem appears to be secondary. As a 
matter of fact the environment quenches the dynamics of these subsystem properties. In a 
classical world this would correspond to a picture of motionless subsystems, like molecules 
in material systems, without fluctuations in the thermodynamic limit. Such a stationarity 
derives from the fact that the expectation value is not a directly observed physical property, 
but an average of infinite measures of the physical property. Nonetheless, the expectation 
value is the standard tool supplied by quantum mechanics for the description of the time 
evolution of physical properties, tool which displays stationarity in the thermodynamic limit. 
It should then be useful to explore the Bolim theory by looking for alternative tools able to 
capture the fluctuation dynamics of parts, like molecules, of a larger isolated system. 


III. PILOT-WAVE THEORY 

Pilot-wave theory is a formulation of quantum mechanics firstly proposed by de Brogli^® 9 
and afterwards rediscovered and fully developed by David Bohm in 195^ 2 3 . It has several 
advantages with respect to more traditional approaches to quantum mechanics, since it leads 
to a full characterization of system constituents as particles having well defined geometrical 
coordinates Q(t) that evolves in time in a deterministic way. Furthermore, it allows a 
description of measurement processes without the need of the wave function’s collapse. 
Beyond these evident benefits, we think that pilot-wave theory is important also because it 
allows a representation of subsystem dynamics which overcomes the stationarity found for 
expectation values in the thermodynamic limit . 

Given the polar representation of the wave function, T( q,t ) = R(q, t)e lS ( q ' t ' > / h , from the 
Schrodinger equation one can derive the time evolution equations for the amplitude R(q , t ) 
and for the quantum phase S(q 1 t). The equivalence between the latter equation and the 
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classical Hamilton-Jacobi equation is the starting point of Bohm’s analysis, which leads to a 
rate equation for particle positions Q(t) controlled by the phase S(q,t) to be interpreted as 
a time dependent held. On this basis, Bohm formulated his theory that can be summarized 
by the following three main assumptions. 

The first concerns the wave function T (g, t ) representing a held which pilots the system’s 
constituents (particles), and evolving in time according to standard Schrodinger equation 
independently of particle positions 

(28) 

Secondly, the instantaneous configuration of the system is specified through the set Q{t ) = 
{Qk(t)}k=i,...,n °f coordinates of all its particles. Like in classical mechanics they evolve in 
time along a trajectory satisfying the following rate equation in the case of particle with 
mass rrik 

d Qk(t) dS(q,t) 

m k -77 = - 

dt d( lk q=Q( t) 

Thus, once the wave function T (g, t) is provided at all times and if the initial configuration 
Q(0) is known, in principle one can derive the trajectory Q{t) for the evolution of the system 
configuration. A physical observable is represented simply by a function A(Q) of the system 
configuration and its value A(Q(t)) at time t is evaluated from the corresponding coordinates 
Q{t). The quantum nature of system dynamics arises from the pilot role of the wave function 
on the coordinate evolution, generating trajectories in general different from those predicted 
by classical mechanics. 

The third assumption concerns the correspondence with the predictions of ordinary quan- 
turn mechanics as specified by the probability density p(g, t ) = | (g, t) | on the coordinates. 
This calls for a statistical description of the system configurations by imagining a swarm of 
trajectories generated by different initial conditions Q(0), but with the same wave function 
as pilot agent. Let us introduce the probability density p 0 (g) for the initial coordinates Q(0). 
Then we can define the probability density p(g, t) of system coordinates at a generic time t, 



P(q,t) ■= J d QoPo{qo)S(q~ Q{q 0 ]t)), (30) 

on the basis of trajectories Q(qo ; t ) starting at Q(qo ; 0) = go, always with the same pilot wave 
function. In this way, a meaningful comparison can be made between p(g, t ) and p(g, t) since 
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both are probability densities on the same variables. The identity of the initial coordinate 
distributions and of the initial quantum probability distribution, 

p(g,o) =p 0 (g) =p(g,o), (31) 


is the third assumption of the Bohm’s theory. This condition, known also quantum equilib¬ 
rium!^, assures the equivalence between coordinates distribution and the squared norm of 
the wave function at all times 2 


p{q,t) =p(q,t) = |tf(g,t)| , 


(32) 


that is, the Born’s rule for the correspondence between wave function and coordinate proba¬ 
bility distribution. In this way, the same predictions are recovered from Bohm’s theory and 
from traditional quantum mechanics. For instance the average of A(Q(t)) on the ensemble 
of trajectories, 

J dq A(q)p(q,t) = A(t), (33) 

becomes equivalent to the expectation value A(t) of Eq. (J2]) for the operator A given as 
function A(q). 


It should be mentioned that distribution Eq. (30) on the ensemble of Bohm’s trajectories 


can be specialized to the subsystem by integration on environment coordinates, 


P S (lS,t)'-= / ilE P(HS,lE,t). 


(34) 


Such a reduced distribution is equivalent to the quantum marginal probability density 


Eq. (24) provided that the quantum equilibrium condition Eq. (32) is satisfied, 


P s (ls, t) = p s (qs, t ) =; p S ’ e, {qs), 


(35) 


with stationarity holding in the absence of fluctuations for large enough systems. 

Despite the original aim of Bohm’s theory to overcome the methodological flaws of a 
traditional quantum formulation, in order to ensure the agreement with predictions of quan¬ 
tum theory, it introduces a new controversial issue in relation to the third assumption and 
the role of Born’s rule. First of all, there are not clear and evident justifications of the 


initial equivalence Eq. (31) between coordinate distribution and squared modulus of the 
wave function. Such a issue has been tackled several times 6 9 10 with different proposals for 
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the mechanisms ensuring the relaxation towards quantum equilibrium Eq. (32) if initially 
Born’s rule is not satisfied. However, we think that the same idea of a swarm of trajec¬ 
tories rises some methodological criticisms. Indeed, by adopting a realistic point of view, 
a particular realization of the quantum system should be described by a single trajectory, 
while the swarm of trajectories should represent an ensemble of realizations of the system 
in correspondence of different initial configurations Q(0). 

We intend in this work to explore the Bohm’s theory at the level of the single trajec¬ 
tory representation of the quantum system. Such an approach is equivalent to the Bohm’s 
theory as long as the correspondence with standard quantum theory through Born’s rule 
is not considered. In this way, two reality elements describe completely the system’s state: 
the coordinates and the wave function held, (Q(t), H/(g, t)). The state evolution is strictly 
deterministic according to equation: 

/ d Q(t) \ 

= x(Q(t)Mq,t)), (36) 


V 


d t 

dv(q,t) 


dt 


with a time independent vector held X(Q, T) derived from Schrodinger equation Eq. (28) 


and from Bolim equation Eq. (29). From a mathematical point of view, the system state is 
represented by an element of the union of the configuration space and of the Hilbert space, 


with a time evolution described by the Cauchy problem Eq. (36) having an unique solution 
for given initial configuration Q 0 = Q( 0) and initial wave function T 0 (g) = T(g, 0). 

In order to recover a probabilistic description from a single Bohm’s trajectory, one has 
necessarily to resort to the statistical sampling of the coordinates during their time evolu¬ 
tion, like in ergodic theory of classical statistical mechanic^. As long as such a sampling 
represents overall effects of system evolution, it is an equilibrium property which should 
in general depend on the constants of motion, that is the populations P determining the 
pilot wave function. The probability density on coordinates q extracted from the sampling 
of a single trajectory Q(t) will be denoted as w eq (q ), keeping implicit the reference to the 
parametric dependence on the population set to deal with a more compact notation. In 
order to perform a meaningful comparison with the distribution obtained from the wave 
function, we shall consider in a multicomponent system the probability density w s,eq (qs) 
on the subsystem S described by coordinates Qs■ As discussed in the previous section in 


relation to Eq. (27), if the isolated system is large enough, the quantum probability density 
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of the subsystem p s (qs, t ) have negligible fluctuations and it can replaced by its equilibrium 
form p s ’ eq (qs), that is a time independent distribution. In other words, by examining a part 
of a much larger system, the quantum distribution can be described by p s,eq (qs), that is a 
time independent function like the coordinate distribution w s,eq (qs ) obtained from a single 
Bohm’s trajectory and, therefore, a meaningful comparison between them can be done. This 
is the objective of the calculations in a model system reported in the next section. 


IV. BOHM’S TRAJECTORIES IN A MULTI-PARTICLE MODEL 
SYSTEM 

In order to compare the single Bohm’s trajectory and the quantum distribution function, 
we have examined the dynamical behavior of a model system of six confined rotors interacting 
through random potentials. The numerical calculations done for a typical situation clearly 
show that a correspondence exists between the Bohm’s coordinate distribution w s,eq (qs ) and 
the equilibrium quantum distribution p s,eq (qg ) for a rotor subsystem. In this section, after 
the presentation of the model system, we discuss the numerical methods employed for the 
calculation of the relevant properties and we illustrate the most relevant results. 


A. The model system 

We shall consider a system of n — 6 identical but distinguishable particles with mass 
m that move on a ring of constant radius R. Such a system is equivalent to n planar 
rotors described by the set of angles q = {qi, q - 2 ,..., q n }, each of them having an inertia 
momentum / = mR 2 . A physical realization of the system could be an ensemble of methyl 
groups rotating about their C — CH 3 bonds. The Hilbert space R, for the i-th rotor is the 
ensemble of periodic functions of the angular coordinate q % , whose Fourier representations 
can be generated by means of the following orthonormal basis set 

e *m 

Xj{ q i) = -j=, (37) 

with integer values for j index. The tensor product of the Hilbert spaces of each rotor, 
% — R\ ® 'H -2 ®. • ■ <8> Rn-, identifies the Hilbert space for the overall system. Such a model 
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system will be described by means of the following Hamiltonian 

n 

H = H m + VO = Y Ht ] + V {r \ 


( 38 ) 


i —1 


where H is the single particle Hamiltonian, while V^ is an interaction potential of random 
type. For the single particle Hamiltonian we use the model of a planar rotor confined by a 
cosine potential with minimum at g,; = n: 

k 2 d 2 


W (0) = 


u 


2 Idqf ' 2 


+ ^(i+COSft), 


(39) 


the parameter u representing the energy barrier at q\ = 0. In the following the parameter 
h 2 /21 will be employed as the energy unit. We intend to analyze the quantum dynamics 
of the system in conditions of significant confinement of the rotors, and to this purpose we 
have selected the potential barrier as u — 300(h 2 /2J). 

The contribution V ^ of the system Hamiltonian has the purpose of producing a dynam¬ 
ical coupling between rotors by means of random interactions typical of material systems. 
Moreover, it assures the rational independence of the Hamiltonian H eigenvalues, property 
which does not hold in the presence of identical single particle Hamiltonians only. The ran¬ 
dom potential has been parameterized as single particle contributions and interaction terms 
between pairs of rotors: 


vf’fe) + \ Y f 1 - Mhfta (40) 

i= 1 i,j=l 

Let us denote with V(0) the periodic function representative of a single particle contribution, 
that is Vy r ' ) (g*) for g* = 6 , or of a two-particle interaction, that is (q, — qj) for g, — qj = 6. 
By means of a gaussian random variable with null average and a given variance cry, a random 
profile is easily generated for its discretized values 14 := V{6k) at (2 L + 1) equally spaced 
angles 6 k = 2nk/(2L + 1) for k = 0,1,... ,2 L. Standard algorithms can be employed to 
produce these random values of the function with statistical properties 

v k = 0, V 2 = a 2 Vr (41) 


where the average is referred to different realizations of the same coefficient. In order to 
recover a continuous function V ( 6) from these random coefficients, we resort to a truncated 
Fourier decomposition 

L 


V(6) = ^ 

l=—L 


(42) 
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with its (2L + 1) components evaluated at the discretized angles 


2 L 


V t = 


2L + 1 


J2 Vke 




( 43 ) 


k =0 


Since an additive constant in the potential does not modify the quantum dynamical prop¬ 
erties, a null value is attributed to the Fourier component Vo, this being equivalent to the 
constraint of a null angular average for functions V(9). 

In conclusion, the previous procedure allows the generation of these random angular 


functions for each contribution of Eq. (40) on the basis of two parameters: the variance ay 


and the number (2 L + 1) of discretized angles. The variance ay controls the strength of the 


random contribution V^ r ’ with respect to single particle Hamiltonians in Eq. (38). In the 


following calculations we shall use an unitary value of this variance in the adopted energy 
units, that is ay = h 2 /21. This corresponds to random potentials with a strength much 
smaller than the confining potential with a barrier height u = 300(ft 2 /2J). In this way, the 
random potential contribution has nearly a perturbation effect so that eigenfunctions and 


eigenvalues of the full Hamiltonian Eq. (38) preserve the main features deriving from single 


rotor contributions. The other parameter L controls the size of the angular correlations in 
the potential, since it determines the distance between two adjacent discretized angles with 
uncorrelated values of the potential. In the calculations we shall use the value L = 100 
because it produces an highly random potential. An angular dependence resembling that of 
a noisy signal is evident form Fig. [l] which displays the potential deriving from a particular 
realization of the (2 L + 1) coefficents 14 for L = 100. For each of the contributions of 


Eq. (40) an independent realization of the random potential V{9 ) is employed. 


B. Numerical methods 


Before to discuss the results for the quantum dynamics of the model system, we summarize 
in this section the numerical procedures employed in the calculations. They concern four 
main issues: i) the solutions for the single rotor Hamiltonian, ii) the eigenvalues and the 
stationary states for the system of interacting rotors, iii) the initial conditions and the time 
dependent wave-function and related properties, iv) the Bohm’s trajectory. 


The eigenfunctions of the single rotor Hamiltonian Eq. (39) are required because their 


tensorial products represent the most convenient basis for the numerical solution of the time 
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Figure 1. Random potential characterized by the parameters ay = fi 2 / 21 and L = 100. 


independent Schroedinger Eq. (JT]) as long as the random potential V ^ is weak. Let us 
denote the eigenvalue problem for the single rotor as 

(44) 

with eigenvalues ordered from below, e m < e m+ll for m = 0,1, 2,.... Of course all the 
rotors have the same eigenvalues and eigenfunctions with the same functional form since 


they share the same Hamiltonian Eq. (39). In order to obtain the relevant eigenfunctions 


and eigenvalues with the required precision, we have generated the matrix representation of 


Eq. (39) on the basis Eq. (37) for \j\ < 20, and diagonalized it by employing the software 
routine Armadillo, a C++ linear algebra library 4 ^. In Table [I] we have reported the lower 
energy eigenvalues, and in Figure [2] the profiles of the corresponding squared eigenfunctions 
\<fm(Qi)\ 2 with the eigenvalues as offset together with the confining potential. In Table |T] we 
have also included the harmonic oscillator eigenvalues resulting from the parabolic approxi¬ 
mation w(l + cos qt )/2 ~ u(qi — 7t) 2 /4 of the rotor potential, in order to attest the differences 
with respect to purely harmonic quantum dynamics. Indeed for increasing levels the dif¬ 
ference between the two eigen-energies clearly emerges. It should be mentioned that the 
numerical diagonalization of the single rotor Hamiltonian supplies not only the eigenvalues 
e m , but also the eigenvectors, that is the coefficients for the expansion of the eigenfunctions 
Vm{qi) on the basis of Eq. ( |37| ). When, in the following, operations on single rotor eigen¬ 
functions <p m (qi) are invoked, implicitly we refer to operations on these linear combinations 
which can be easily encoded in computer programs. 


Given the numerical solutions of Eq. (44), one can employ the following basis for the 
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Figure 2. Squared modulus of eigenfunctions c p m (qi ) of the 7 /j 0 - 1 Hamiltonian with u = 300(h 2 /2 1) 
as the intensity of the conhning potential u( 1 + cos(< 7 i ))/2 


Table I. Low energy eigenvalues e m of single rotor Hamiltonian for the potential barrier u = 
300(h 2 /2J). The corresponding harmonic oscillator eigenvalues are reported between parenthesis. 


m 

e m /(h 2 /2J) 

0 

8.597 (8.660) 

1 

25.664 (25.981) 

2 

42.472 (43.301) 

3 

59.015 (60.622) 

4 

75.286 (77.942) 

5 

91.278 (95.263) 

6 

106.982 (112.583) 

7 

122.390 (129.904) 

8 

137.491 (147.224) 

9 

152.275 (164.545) 


Hilbert space TL of the overall system 

n 

10 = 01 Wi(9t)>, (45) 

i=l 

where l := (Zi, Z 2 , ••• ,/ n ), and each index identifies the eigenfunction of the cor¬ 

responding i-t\i rotor with m = l t . The basis elements |Z) are eigenfunctions of the model 
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Figure 3. Numerical results of energy eigenstates for two truncation parameters E^) = 
154/(h 2 /2I) (blue crosses) and E^) = 171/(Ti 2 /2I) (red crosses). The states with comparable 
energy are labeled according to the polyad quantum number, from P=0 (ground state) to P=7. 
The last part of polyads P = 5 and P = 6 are magnified in the insets in order to show the effects 
of the truncation parameter. 


system Hamiltonian in the absence of the random potential 

n 

ff<°> io = £f ) |/>, £f> = y>„ 


(46) 


i=l 


and they are conveniently ordered according to the corresponding energies e[°\ As long as 
the random potential V ^ acts like a perturbation, the diagonalization of the full Hamil¬ 
tonian is influenced mainly by the coupling between basis elements with nearby values of 
E\ 0 \ and this allows an efficient truncation of the Hamiltonian matrix representation. In 
practice, one considers all the basis elements with Ef 0 '* less than a given truncation energy 
cutoff E^). Then the matrix representation of the full Hamiltonian is generated in order to 
perform the diagonalization by means of the software Armadillo. 

The organization of these energy levels in well separated multiplets is evident in analogy 
to the polyads describing molecular vibrations (see 41 42 and references therein). Since in the 
harmonic approximation the oscillators for the confined rotors are degenerate, the polyad 


quantum number classifying the basis elements Eq. (45) is given as P = Y^=i h with values 
P = 0 (ground state), P = 1 (6 states), P = 2 (21 states), and so on. As long the random 
potential V ^ is weak, the corresponding Hamiltonian eigenfunctions \Ek) are substantially 
reproduced by linear combinations of basis elements with a given polyad quantum number 
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P, with only perturbational contributions from the other polyads. Therefore the polyad 
quantum number can be used to classify also the eigenvalue multiplcts, as done in Fig. [3] 


The comparison in Fig. [3] of the numerical eigenvalues obtained with two values of pa¬ 
rameter E^} allows one to evaluate the effects of matrix truncation. Notice that the chosen 
values of E$ leads to a complete inclusion of the selected polyads in the truncated matrix 
representation. The results with matrix representation for E^/ = 15A/(h 2 /2I) (N = 924), 
polyads from P = 0 to P = 6) will be employed as the reference for the calculation of time 
dependent properties of the model system. Their accuracy has been checked by comparison 
with the larger matrix obtained for = 171/ (h 2 /2I) which includes a further polyad. 
Such a matrix enlargement, besides introducing new eigenvalues (i.e., the P = 7 polyad), 
produces a change of about 0, 04% for the upper energy eigenvalues, and smaller variations 
for decreasing energy (see the insets of Fig. [3]). Such a behavior agrees with the pertur¬ 
bational contribution by the random potential surely it has strong effects within a 

polyad in the presence of degenerate or nearly degenerate zero-order energies e[°\ but it 
has weak effects on the coupling of states belonging to different polyads with well separated 
values of e[° . These informations allows us to conclude that with the truncation parameter 
E$ = 154(7 2 /2J) (N = 924) we get numerical results with errors at most of 0,04%. As 
a matter of fact the accuracy of the data employed in the calculation of time dependent 
properties is much better, as long as we shall use an active space including up to P = 5 
polyad whose eigenvalues deviate from those obtained with the larger matrix by 0, 004% at 
most. The final results of this computational task is the ensemble of eigenvalues E *. and 


eigenfunctions 17%), the latter specified as linear combinations of basis elements Eq. (45) 
through coefficients (l\Ek)), for the time independent Schroedinger Eq. ([7]). 


Once the eigenstates and the energy eigenvalues are obtained and the active space is 
identified on the basis of the cutoff energy E max , the time dependent wave-function has 
to be evaluated. Thus the initial quantum state has to be chosen according to the set 
of populations P and the set of initial phases a(0) within the active space. Since the 
phases are homogeneously distributed®, they are simply selected at random within their 
domain. Also for the populations a random choice is performed but, in order to preserve their 
normalization, by means of suitable set of auxiliary parameters homogeneously distributed 
in the (0,1] domain according to procedure discussed in Ref. 30 43 ®. Given these initial 
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conditions, the wave-function at an arbitrary time is specified as 


N 




( 47 ) 


k =1 


where N is the dimension of the active space. For the calculation of the reduced density 
matrix, reference is made to the first rotor, qs = Qi, so that its matrix elements on the basis 


of single rotor eigenfunctions Eq. (44) can be specified as 


&m,m' {t') ■ i}Pm \) | Pm ') 

=E ik-d mm <*wio = 


1,1' \i =2 


(48) 


=ee n | m^k) (Ek'\l ) X 

k,k' 1,1' \i=2 / 

x v /p fc^'e _l[Qfe(0)_Qfe,(0)+(Efe_E '=' )f/?i] . 

The same equation with the constraint k — k' in the summations on the r.li.s. can be 
employed to evaluate the elements of the equilibrium density matrix. Given the 

reduced density matrix, also the marginal quantum distribution of the subsystem (the first 
rotor) is recovered according to Eq. 


P S (QS,t ) = (49) 

771,77?/ 

and the equilibrium distribution as well by inserting the equilibrium density matrix elements 


P {.Qs) ^ ^ ®m,m' ( pm{,qs') t Pm'(SLs) • 

771, m' 


(50) 


By specifying the eigenfunctions | Ek) in Eq. (47) as linear combinations of the basis functions 


Eq. (45), one gets for a given time the explicit dependence on the coordinates of the wave 


function, 4/(g, t), and of both the amplitude R(q,t) and the phase S(q,t) as well. 

For the computation of the trajectory of the rotors, we adopted the Runge-Kutta 
method 45 at the 4 th order to solve the Bohm equation of motion Eq. ( |29| ). We employed a 
time step A t = 0.01(47r//ft) that assures a good approximation to the calculated trajectory 
from the point of view of its statistical properties. In particular we have evaluated the 
correlation function G(t) of the rotor angle Qs 


G(t) := AQ s (t)AQ s (t + r), 


(51) 
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with A Qs{t) = Qs(t) — Qs(t), that we calculate from the discretized time average along the 
trajectory: 

1 M 

G(r) -jjY. &QsUAt)AQ s QAt + r), (52) 

V1 3=0 

where M is the number of sampling points that depends on the length of the examined 
trajectory. 

Finally the distribution w s,eq (qs ) of the rotor coordinate along its trajectory has to be 
evaluated. In practice we have calculated its discretized counterpart by dividing the domain 
0 < Qs < 2vr of the rotor angle into 10 4 equally spaced intervals. The probability density is 
recovered from the fraction of time spent by the rotor in each interval during its evolution. 
In order to check that the length of the trajectory is sufficient, we have verified that the 
resulting distribution is not modified by a further evolution. 

One might wonder whether the numerical procedure for the calculation of the Hamilto¬ 
nian eigenvalues and eigenfunctions, which provides always approximate results, affects the 
behavior of the computed trajectory. If this is the case, then in the comparison between 
the quantum distribution p s,eq (qs) for the subsystem and the distribution w s,eq (qs) on the 
Bohm’s coordinate, one should consider explicitly the influence of the errors introduced 
by the numerical diagonalization. Let us denote with E[ app ^ and \E[ app ' > ) the approximate 
eigenvalues and eigenfunctions computed numerically. The computed wave-function derives 
from the linear combinations of these approximate eigenfunctions, and it is a solution of the 
Schroedinger equation for the Hamiltonian 


N 

ff(app) ._ | E^ app> (E k 

k =1 


(53) 


instead of the assumed model Hamiltonian Eq. (38). Correspondingly the phase function 


S(q,t ) and the resulting Bohm’s trajectory is exact for the quantum problem described by 
the Hamiltonian H^ a PP\ In conclusion, the unavoidable errors introduced by the numerical 
diagonalization are formally equivalent to a slight modification of the system Hamiltonian. 
Of course the Bohm’s trajectory is also affected by the numerical errors in the integration of 


the differential equation (29), but their effects can be easily controlled by checking that the 


coordinate distribution iv s,eq (qs), and the coordinate correlation function Eq. (51) as well, 
does not change by decreasing the integration time step. 
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C. Dynamical properties 


The selected model system and its Hamiltonian with u = 300(fr 2 /2J) as the barrier of 
the confining potential, is compatible with different thermal states depending on the cut-off 
energy E max of the RPSE. We have selected the value E max = 139(fr 2 /2J) which corresponds 
to an active space of dimension N = 462 including the polyads from P = 0 to P = 5 and 
exluding the other polyads (see Fig. [3]). With such a choice we deal with a state having a 
significant distribution between the ground state and the excited states of the single rotor, 
as witnessed by the subsystem reduced density matrix represented on the basis of single 


rotor eigenfunctions ip m (qs) of Eq. (44). The equilibrium reduced density matrix calculated 
according to the methods illustrated in the previous section is nearly diagonal, and the 
diagonal components are reported in Table [TTj The calculated off-diagonal elements cf m ,m' 
are less than 1/1000 in magnitude with respect the associated diagonal elements and 


& m' ,m' • 


The decrease of the diagonal elements cr m ^ m with the single rotor energy e Tl 


see 


Table [n| might suggest a canonical form oc exp(— /3e m ) but this is not the case. In 

order to provide evidences about it, we have derived the hypothetical canonical thermal 
coefficient (5 = 0.0376/(fr 2 /2J) from the ratio aqq/bAo and then the corresponding elements 
of the canonical density matrix, which are reported between parenthesis in Table |TTJ The 
deviations with respect to the numerical values of cf rn . m clearly emerge, particularly for the 
upper energy elements, and this points out that the size of our model system (six interacting 
rotors) is not large enough to ensure the thermodynamic limit. On the other hand the rather 
small differences for the more populated states (say, for m = 0,1,2) allows us to drawn 
the conclusion that the state of the system resembles that of thermodynamic equilibrium. 
Furthermore, the significant mixing of single rotor eigenstates clearly emerges from the data 
of Table [TTJ since excited states (i.e., the elements with m/0) contribute by nearly 50% to 
the reduced density matrix. 


As explained in the previous section, the instantaneous reduced density matrix allows 
the calculation of the time dependent quantum distribution p s (qs,t ) of the subsystem (the 
first rotor). The profiles of such a distribution are reported in Fig. [I] for a selected sample of 
times. As the reference for the visualization of its change with the time, in the same Figure 
we have plotted also the equilibrium quantum distribution p s ’ eq (qs) calculated according 
to equilibrium reduced density matrix In Section II we have shown that in the 
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Table II. Diagonal elements of the equilibrium reduced density matrix, with their canonical values 
reported between parenthesis. 


m 

& m,m 

0 

0.536 (0.475) 

1 

0.282 (0.250) 

2 

0.127 (0.133) 

3 

0.0431 (0.0712) 

4 

0.0122 (0.0386) 

5 

5.15 10“ 4 (0.0211) 

6 

3.61 10“ 7 (0.0117) 



Figure 4. Equilibrium marginal density distribution p s,eq (qs ) (black thick line) and marginal den¬ 
sity distributions p s (qs,t) (colored thin lines) at some selected times. The marginal distributions 
are referred to the first of the 6 rotors in our model system. 


thermodynamic limit, i.e., when the number of interacting components is large enough, the 
fluctuations of p s (qs,t ) become negligible and then p 3,eq (q s ) would reproduce the quantum 
distribution function at all times. The data in Fig. [4] clearly show that this is not the case 
in our model system as long as time dependent deviations from p s,eq (qs ) are evident. On 
the other hand these deviations have a comparably low magnitude, so that the equilibrium 
distribution p s ' eq (qs) can be considered as representative, at least approximately, of the 
instantaneous quantum distribution p 3 (qs,t). 
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Figure 5. Time evolution of the Bohm coordinates (drawn with different colors) of the six rotors 
of the model system. 



Figure 6. Correlation function G(r) of the first rotor coordinate. 


Having characterized the main quantum properties of the subsystem, we examine now the 
behavior of the Bohm’s coordinates. By employing the procedure illustrated in the previous 
section, the trajectories of the angular coordinates Qi(t) of the six rotors have been computed 


according to Eq. (29) by choosing Qj(0) = n as initial conditions in correspondence of the 
bottom of the rotor confining potential. In Fig. [5] we have represented with different colors 
the trajectories of all the rotors within the time window 0 < t/(4irl/h) < 5. Each rotor 
coordinate follows a strongly confined dynamics with limited excursions about the potential 
minimum. The time evolution of each rotor coordinate seems that of a fluctuating signal 
loosing correlation with time, somehow like in the brownian motion. To verify this feature, 


we have computed the correlation function G(r) Eq. (51) which is displayed in Fig. kSl As 
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Figure 7. Equilibrium marginal quantum distribution p s,eq (qs ) and marginal density distribution 
w s,eq (qs ) of the Bohm coordinate for the first rotor. 

expected on the basis of the behavior of the trajectories, the correlation vanishes with a 
rather short correlation time of order t c /(4itI/ h) ~ 4 supporting the analogy with brownian 
motion. It appears that the phase function S(q,t) due to the wave-function generates a 
fluctuating evolution of the Bohm’s coordinates, which leads to a fast loss of correlation. 

If the Bohm’s coordinate Qs of the subsystem is considered as a stochastic process, 


than its properties are naturally characterized by the correlation function Eq. (51) and 


its equilibrium distribution w s,eq (qs)- The expected confinement of the rotor angle clearly 
emerges from such a distribution which is displayed in Fig. [TJ In the same Figure we have 
plotted also the equilibrium quantum distribution for the sake of comparison. The two 
equilibrium distributions, the one p s,eq (qs) deriving from the evolution of the wave-function, 
and the other w s,eq (qs ) calculated from a single Bohm’s trajectory, result to be very close 

w S ’ eq (q s ) ~p S ’ eq (q s ). (54) 

Notice that the loss of correlation along the trajectory implies that the distribution w s,eq (qs) 
is independent of the choice of the initial values Q(0) of the Bohm’s coordinates. 


It should be stressed that the correspondence Eq. (54) cannot be considered as a general 


property for all quantum systems. Indeed one can use a single rotor system as a counterex¬ 


ample where Eq. (54) does not hold. If the same previous procedure is applied to an isolated 


conEned rotor, with the same potential of our model system, by choosing an active space of 
dimension IV = 2 in order do deal with a wave-function with a nearly 50% probability of the 
ground state like for the reduced density matrix of Table |TT[ one obtains two very different 
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Figure 8. Quantum equilibrium distribution p eq (q ) and the distribution w eq {q ) of the Bohm coor¬ 
dinate for the model system with a single rotor. 


equilibrium distributions like those displayed in Fig. [8] It should be mentioned that the 
asymmetry on the profile of w eq (q ) derives from the difference of the randomly chosen pop¬ 


ulations of the two quantum states. We conjecture that the correspondence Eq. (54) found 
in our model system is a consequence of the multi-particle interactions which are absent in 
the single rotor system. 


Finally we emphasize that the correspondence (54) is not an accidental result of particular 
conditions employed for the calculation in our model system. As a matter of fact we have 
similar evidences of the correspondence from calculations in other conditions, for instance 
by using a confining potential with lower strength or even in the absence of the confining 
potential. 


V. BOHM COORDINATES AS MARKOV STOCHASTIC VARIABLES 


As reported in the previous section, explicit calculations with a many body model system 
suggest that, even by considering pilot wave theory at the level of single Bohm’s trajectory, 


a correspondence exists according to Eq. (54) between probability density determined by 


the wave-function and the coordinate distribution derived from the trajectory. This is an 
important result since it allows a connection between the standard quantum theory and 
the geometrical description of system evolution through a trajectory, without the need of 
the Born’s rule and the corresponding swarm of trajectories. On the other hand, such 
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a connection has a methodological role different from the Born’s rule. In particular we 
emphasize the following three features. 1) It is a correspondence concerning the subsystem 


only, while the Born’s rule Eq. (32) deals with the overall isolated system. 2) Its validity has 


to be restricted to the case of negligible fluctuations on the quantum distribution Eq. (23) 


which then can be replaced by the equilibrium quantum distribution Eq. (24). Only in this 


case the quantum distribution becomes time independent and, therefore, it can be compared 
with Bhorn’s coordinate distribution which, by definition, is time independent. 3) Such a 
correspondence in general is not exact, since there are evident counterexamples (Figure [8j) . 
The calculation results simply suggest that quantum distribution and Bohm’s trajectory 
coordinate distribution are close in suitable conditions. On the contrary the Born’s rule 


Eq. (32) is exactly veriEed once the third Bohm’s assumption Eq. (31) is introduced. 


Besides these considerations, an important issue naturally arises: can the correspondence 


Eq. (54) End a support beyond the evidences resulting from calculations with specific model 


systems? In other words, can Eq. (54) be derived under particular conditions? A positive 


answer is found if the evolution of subsystem coordinate in the single Bohm’s trajectory fol¬ 
lows a stationary Markov process for a stochastic variable®. A stationary Markov process is 
completely characterized by the equilibrium distribution w s,eq (qs ) and the conditional prob¬ 
ability distribution w s (qsfi\qs, T )- The former is obtained from the sampling of subsystem 
coordinate Qs(t) along a single Bohm’s trajectory, as we have done in our model system. 
The latter requires the sampling of the correlation of coordinates Qs(t ) and Qs(t + r ) at two 
times separated by r, and it should satisfy the constraint of correlation loss at long enough 
times: 

lim w s (q Sfl \q s ,T) = w s ’ eq (q s ). (55) 

T—>-+00 

The distributions characterizing the Markov process observed in a single trajectory, can 
be used to describe also the probability density arising from an ensemble of trajectories. Let 
us denote with p s (qs,t ) the probability density on the coordinate for such an ensemble of 
trajectories. Given the initial distribution p s (qs, 0), the probability density at any time can 
be evaluated on the basis of the correlation function w s (qsfi\<ls, t), 


P S (ls,t)= / il'Js-ji P S (is,o,0)w s (q s ,o q s ,t), 


(56) 


and, according to Eq. (55), it will relax to the equilibrium distribution w s,eq (qs) at long 
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enough times 


lim p s (q s ,t) = w s ’ eq (q s ). 

T—>■+OO 


( 57 ) 


Let us now recognize the conditions under which the distribution on the trajectories is 
stationary, that is p s (qs,t) is time independent. Stationarity means that the equivalence in 


Eq. (57) must be verified at all times, 


p a (qs,t)=w s ’”(qs), 


and, therefore, the initial distribution p s (qs, 0) = w s,eg (qs ) is the unique condition leading 


(58) 


in Eq. (56) to a time independent distribution. 


Let us apply these results to the ensemble of trajectories generated according to the third 
Bohm’s assumption, and with the further conditions: a) stationary quantum distribution for 
the subsystem, p s (qs,t ) = p S ' eq (qs), b) stationary Markov process for subsystem coordinate 
in a Bohm’s trajectory. Then the subsystem probability density, according to the Born’s 


rule Eq. (32), is equivalent to the stationary quantum distribution 


p a (qs,t) = p a (qs,t) = p s ' e, (qs), 


(59) 


because of condition a). On the other hand, because of condition b), the same probabil¬ 
ity density can be computed according to Eq. (56) with p s (qs, 0) = p s,eq (qs) as the initial 


distribution. However, such a probability density has only one stationary form specified 


by Eq. (58). One can conclude that, as long as the two conditions a) and b) are satisfied 


and, therefore, Eqs. (58) and (59) are holding simultaneously, for the subsystem the quan¬ 


tum equilibrium distribution and the coordinate distribution in a Bohm’s trajectory are 
equivalent, 

p s,eq {qs ) = w s,eq (qs)- (60) 

This is the important result of the previous analysis which, however, is conditioned by the 
validity of assumptions a) and b). In section II we have analyzed the fluctuations of the 
quantum distribution for the subsystem, by showing that they vanish in the limit of an 
infinite size environment. This points out that in finite but large enough systems, condition 
a) is satisfied only approximately, and the same type of validity should be attributed to 


the equivalence Eq. (60). At this stage a specific analysis about the general validity of the 


description of subsystem coordinate as a Markov process is lacking even if, in analogy to 
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classical brownian motion, one might conjecture that such a feature is determined by the cou¬ 
pling amongst several degrees of freedom. On the other hand the model results reported in 
the previous section suggest that for systems characterized by random interactions amongst 


its components, the subsystem evolution leads to distributions approximating Eq. (60). 


VI. CONCLUSION 

We have considered the system of six confined rotors as a model to test the representation 
of quantum systems by the single Bohm’s trajectory. For the subsystem identified with one 
rotor, the others playing the role of the environment, we have found the following main 
results from the numerical solution of the pilot wave theory: 1) the marginal quantum 
distribution derived from the wave-function is nearly stationary, 2) the Bohm’s coordinate 
evolves like a randomly fluctuating signal with a clear loss of correlation with the time, 
3) the rather close correspondence between the marginal quantum distribution and the 
distribution of the Bohm’s coordinate. We stress the interest of the last result in relation to 
the methodological status of the pilot wave theory. If a correspondence exists between the 
quantum distribution derived from the wave-function and the distribution of the Bohm’s 
coordinate along a trajectory, albeit at the level of the subsystem, then the ensemble of 
trajectories together with the Born’s rule for their initial distribution is not mandatory to 
establish a connection between standard quantum theory and the particle’s configuration 
of the pilot-wave theory. A more direct picture of material systems would then be derived 
on the basis of a single realization of both the quantum state (the wave-function) and the 
particle’s configuration (the Bohm’s trajectory). 

On the other hand we emphasize that such a correspondence is presently a conjecture 
even if supported by the numerical results for a particular model system. The existence 
of conditions assuring the validity of the conjecture remains still an open issue. We were 
able to verify it only under the hypothesis that the Bohm’s coordinate behaves as a Markov 
stochastic process, as shown in the previous section. 

Finally we would like to comment on the implications of the near stationarity of the 
marginal quantum distribution for the subsystem, as shown in Fig. |4j This is strictly a 
consequence of the vanishing of fluctuations of the reduced density matrix in the thermody¬ 
namic limit as analyzed in Ref. 25 . A direct relation exists also with the typicality analyzed in 
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Ref. 381 even if in those theories the effects of fluctuations in time are not separated from the 
distribution within the ensemble. At any rate, a static picture of the subsystem properties 
would be implied from these results, at odds with the opposite image of an ever fluctuat¬ 
ing world as suggested by the classical mechanics. We think that the pilot wave theory, in 
the single trajectory approach, leads to a solution of these contradictory representations. 
Indeed, as in displayed in Fig. [5j the fluctuating evolution of the Bohm’s coordinate, very 
much like for a confined Brownian particle, results compatible with the nearly stationarity 
of the marginal quantum distribution. 
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